Spin-valley interplay in two-dimensional disordered electron liquid 
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We report the detailed study of the influence of the spin and valley splittings on such physical 
observables of the two-dimensional disordered electron liquid as resistivity, spin and valley suscepti- 
bilities. We explain qualitatively the nonmonotonic dependence of the resistivity with temperature 
in the presence of a parallel magnetic field. In the presence of either the spin splitting or the valley 
splitting we predict novel, with two maximum points, temperature dependence of the resistivity. 
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I. INTRODUCTION 

Disordered two-dimensional (2D) electron systems 
have been in the focus of experimental and theoretical 
research for several decades. 1 Recently, the interest to 2D 
electron systems has been renewed because of the exper- 
imental discovery of metal-insulator transition (MIT) in 
a high mobility silicon metal-oxide-semiconductor field- 
effect transistor (Si-MOSFET)^ Although, during last 
decade the behavior of resistivity similar to that of 
Ref. 011] has been found experimentally in a wide va- 
riety of two-dimensional electron systems^ the MIT in 
2D calls still for the theoretical explanation. 

Very likely, the most promising framework is provided 
by the microscopic theory, initially developed by Finkcl- 
stein, that combines disorder and strong electron-electron 
interaction^ Punnoose and Finkelstein^ have shown pos- 
sibility for the MIT existence in the special model of 2D 
electron system with the infinite number of the spin and 
valley degrees of freedom. The current theoretical re- 
sults^ do not support the MIT existence for electrons 
without the spin and valley degrees of freedom. There- 
fore, it is natural to assume that the spin and valley de- 
grees of freedom play a crucial role for the MIT in the 
2D disordered electron systems. 

Usually, in the MIT vicinity, from the metallic side, 
i.e., for an electron density higher than the critical one, 
and at low temperatures T <C r t 7 the initial increase 
of the resistivity (ft) with lowering temperature is re- 
placed by the decrease of p as T becomes lower than 
some sample specific temperature.— Here, r tr denotes the 
elastic scattering time. This nonmonotonic behavior of 
the resistivity has been predicted from the renormaliza- 
tion group (RG) analysis of the interplay between dis- 
order and electron-electron interaction in the 2D disor- 
dered electron systems As a weak magnetic field B 
is applied parallel to the 2D plane, decrease of the re- 
sistivity is stopped at some temperature and p increases 
again^ Further increase of B leads to the monotonic 
growth of the resistivity as temperature is lowered, i.e., 
to an insulating behavior, in the whole T-range. These 
experimental results suggest the significance of the elec- 
tron spin for the existence of the metallic phase in the 
2D disordered electron systems. 



As is well known, in both Si-MOSFET- and n-AlAs 
quantum wellii 2D electrons can populate two valleys. 
Therefore, these systems offer the unique opportunity 
for an experimental investigation of an interplay between 
the spin and valley degrees of freedom. Recently, using a 
symmetry breaking strain to tune the valley occupation 
of the 2D electron system in the n-AlAs quantum well, as 
well as a parallel magnetic field to adjust the spin polar- 
ization, the spin - valley interplay has been experimen- 
tally studiedj 12 ' 13 However, the electron concentrations 
in the experiment were at least three times larger than 
the critical onepii Therefore, the spin - valley interplay 
has been studied in the region of a good metal very far 
from the metal-insulator transition. 

In the present paper we report the detailed theoreti- 
cal results on the T-behavior of the 2D electron system 
with two valleys in the MIT vicinity. In particular, we 
study the effect of a parallel magnetic field and/or a val- 
ley splitting (A„) on the transport, and the spin and 
valley susceptibilities. We find that in the presence of 
either the magnetic field or the valley splitting the metal- 
lic behavior of the resistivity survives down to the zero 
temperature^ For example, this result implies that at 
B = the metallic p{T) dependence can be observed ex- 
perimentally at temperatures T >C A„. Only if both the 
magnetic field and the valley splitting are present, then 
the metallic behavior of the resistivity crosses over to the 
insulating one. Next, we predict novel, with two maxi- 
mum points, T-behavior of the resistivity in the presence 
of the magnetic field and/or the valley splitting. Finally, 
we find that as T vanishes the ratio of the valley sus- 
ceptibility (xv) to the spin one (x s ) becomes sensitive 
to the ratio of the valley splitting to the spin one. At 
high temperatures the ratio Xv/Xs is temperature inde- 
pendent and can be chosen equal unity. If the spin split- 
ting is larger (smaller) than the valley splitting, then at 
low temperatures the ratio Xv/Xs < (>)!■ If the spin 
and valley splittings are equal each other, then the ratio 
Xv/Xs = 1 as temperature vanishes. 

The presence of the parallel magnetic field and the 
symmetry-breaking strain introduces new energy scales 
A s = gLfJ-BB and A„ in the problem. Here, gi, and 
Pb stand for the g-factor and the Bohr magneton, re- 
spectively. Let us assume that the following conditions 
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hold: A„ <C A s <c l/r tr . In addition, a magnetic field 
B± > T/(De) is applied perpendicular to the 2D elec- 
tron system in order to suppress the Cooper channel. 
Here, e and D denote the electron charge and diffusion 
coefficient, respectively. Due to the symmetry breaking, 
the spin and valley splittings set the cut-off for a pole 
in the diffusion modes ( "diffusons" ) with opposite spin 
and valley isospin projections. In the temperature range 
A s «T«T te ' this cut-off is irrelevant and the 2D elec- 
tron system behaves as if no symmetry breaking terms 
are applied. The temperature behavior of the resistiv- 
ity is governed by one singlet and 15 triplet diffusive 
modest At low temperatures A„ < T < A s , eight dif- 
fusive modes with opposite spin projections do not con- 
tribute. Then, the p(T) dependence is determined by 
the remaining one singlet and seven triplet modes. As 
we shall demonstrate below the behavior of the resistiv- 
ity can be either metallic or insulating. Surprisingly, we 
found that the seven triplet diffusive modes are not equiv- 
alent. They have to split into two groups of six and one 
modes for the spin susceptibility be T-independent. For 
temperatures T <C A„, next four diffusive modes with 
opposite isospin projections become ineffective. In this 
case, the temperature dependence of the resistivity is de- 
termined by one siglet and three triplet diffusive modes. 
Although, the number of the remaining diffusive modes 
corresponds formally to single- valley electrons with spin, 
the p(T) behavior is insulating. 

The paper is organized as follows. In Section |TT] we 
introduce the nonlinear sigma model that describes the 
disordered interacting electron system. Then, we con- 
sider the short length scales at which the system has 
SU(4) symmetry in the combined spin and valley space 
(Sec. IHIj) . The behavior of the system at the interme- 
diate and long length scales is studied in Sec. [IV] and 
Sec. El respectively. We end the paper with discussions 
of our results and with conclusions (Sec. IVIjl . 



II. FORMALISM 

A. Microscopic Hamiltonian 

To start out, we consider 2D interacting electrons 
with two valleys in the presence of a quenched disor- 
der and a parallel magnetic field at low temperatures 
T < t£ x . We assume that the perpendicular magnetic 
field B± > T/(De) is applied in order to suppress the 
Cooper channel. Using one electron orbital functions, we 
write an electron annihilation operator as 



Vv(R) = £ W(r)<p(z)e 



irzQ/2 



(1) 



where z denotes the coordinate perpendicular to the 2D 
plane, r the in-plane coordinate vector, and R = r + ze z . 
The subscript r enumerates two valleys and tf>° is the 
annihilation operator of an electron with the spin and 
isospin projections equal a/2 and t/2, respectively. Let 



us assume that the wave functions ip(z) exp(HQz/2) 
are normalized and orthogonal with negligible overlap 
/ dz f 2 (z) exp(iQz). The vector Q = (0,0, Q) corre- 
sponds to the shortest distance between the valley min- 
ima in the reciprocal space: Q ~ a^~ a ] , with ai a t being the 
lattice constant 

In the path-integral formulation 2D interacting elec- 
trons in the presence of the random potential V(r) are 
described by the following grand partition function 



with the imaginary time action 



S = 



l/T 



dt 



-W(r,t)d t iPZ(r,t)-H -U dis -H 



int 



(2) 



(3) 



The one-particle Hamiltonian 

V 2 , A s 



Ho = / dr^(r) 



2m„ 



A v 



(4) 



describes a 2D quasiparticle with mass m e in the pres- 
ence of the parallel magnetic field and the valley splitting. 
Here, /i denotes the chemical potential. Next, 



H dis = J dr^>)Wr)^ 2 (r) 
involves matrix elements of the random potential: 



,(r) = J dzV(R)(p 2 (z)e 



i(r 2 -T 1 )Qz/2 



(5) 



(6) 



In general, the matrix elements V TlT2 (r) induce both the 
intravalley and intervalley scattering. We suppose that 
V(R) has the Gaussian distribution, and 

(V(R)> = 0, (V(Ri)V(RiO) = W(|n - r 2 |, \ Zl - z 2 \), 

(7) 

where the function W decays at a typical distance d. If 
d is larger than the effective width of the 2D electron 
system, i.e., d 3> [ / dz ip A (z)]~ l , then one can neglect the 
^-dependence of V(R) under the integral sign in Eq. ©. 
In this case, the intravalley scattering survives only: 

(8) 



In the opposite case, d <C [/ dz<p 4 (z)} 1 , one find; 



,16 



<VWri)VWr 2 )} 

+5. 



W(|ri-r 3 |,0) (9) 



TirAs^n ~T 2 I^(|ri - r 2 |,2Qr 2 ) 



dz if * (z), 



where W(r,Q) = J dzW(r, z) exp(iQz). The other cor- 
relation functions, e.g. with t\ = — r 2 and T3 = T4, vanish 
due to integration over (21 + z 2 )/2 coordinate. It is the 
last term in Eq. ^ that contributes to the intervalley 
scattering rate 1/t v . Assuming Q" 1 <C d, one can ne- 
glect the intervalley scattering rate in comparison with 
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the intravalley scattering rate 1/tj ~ W(r,0). At last, 
allowing for a low electron concentration n e in 2D elec- 
tron systems, we consider the case when the following 
inequality holds, n e d 2 < 1. Then, both Eqs. © and © 
read 

(K^OOKar^)) = -^^ r 5 TlT J T3Ti 8{Tl - r 2 ), (10) 



where S a represents the free electron part 2 ^ 
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Tt(VQY 



(17) 



2ttv I d 2 rdz 1 dz 2 W(r,\z 1 - z 2 \)ip 2 (z 1 )ip 2 (z 2 ), 



Q- 1 <^d,[J i P 4 (z)dz}- 1 < n- 1 / 2 . 



(11) 



Here, v is the thermodynamic density of states. Under 
conditions (jlip . the interaction part of the Hamiltonian 
is invariant under global SU(4) rotations of the electron 
operator in the combined spin-valley space: 



.„ e 2 f , , ^(rQ^^n)^^)^^) 

/tint = — / aridr 2 ; : . 

2e./ ri-r 2 



(12) 

A dielectric constant of a substrate is denoted as e. The 
low energy part of TCi n t can be written as ^ 17 ' 18 i 19 

Hint = g y rfridr 2 p(ri)r a (r x - r 2 )/a(r 2 ) (13) 
+m a (r 1 )r t (r 1 -r 2 )m a (r 2 )] (14) 

p(r)=^^(r)C(r), (15) 

(TT 



where 



Here, cr.^ denotes the mean-field conductivity in units 
e 2 /h. The symbol Tr stands for the trace over replica, 
the Matsubara frequencies, spin and valley indices as well 
as integration over space coordinates. 
The term£ 

S F = 4nTzTrr){Q- k)+-nTT f d 2 r ^ tr/£Qtr/^ n Q 

an 

-ttTT 2 f d 2 r^(tr/«Q)®(tr/f„Q). (18) 

"* an 

involves the electron-electron interaction amplitudes 
which describe the scattering on small (T) and large 
(r 2 ) angles and the quantity z originally introduced by 
Finkelstein 5 which is responsible for the specific heat 
renormalization.-2i The interaction amplitudes are re- 
lated with the standard Fermi liquid parameters a o 5 ' 17 i 19 
T 2 = -zFgftl + F£), 4r = T 2 + zFg/(l + Fg), and 
z = nv* /2 where v* — m*/(2n) with m* being the ef- 
fective mass. The case of the Coulomb interaction corre- 
sponds to the so-called "unitary" limit,— Fq — > oo. 

The symbol tr involves the same operations as in Tr 
except the integration over space coordinates, and tr A ® 
ixB = A^ZB^^. The matrices A, r) and F h are 
given as 



A^ 2 =signK)<WS Q/3 £ ClC2 

( ma/3;CiC2 _ X , A Q TA/ 3 TACif2 

\ 1 k)nm — Vn-m,kU " O 



(19) 



Here, T s (q) = U {q) + Fk° / (Av) involves the long-range In . the absence of A s and A„, the action S„ + S F 



part of the Coulomb interaction U{q) — 2iTe 2 /(qe) and 
T t (q) = Fq/(4v). Quantities Fq and Ffi are the stan- 
dard Fermi liquid interaction parameters in the singlet 
and triplet channels, respectively. The matrices t a with 
a = 1, ... 15 are the non-trivial generators of the SU(A) 
group. 



B. Nonlinear sigma model 

At low temperatures, Tr tr <C 1, the effective quantum 
theory of 2D disordered interacting electrons described 
by the Hamiltonian (J3]) is given in terms of the non-linear 
cr-model. This theory involves unitary matrix field vari- 
ables Qmn-Vrrf 2 ( r ) which obey the nonlinear constraint 
Q 2 (r) = 1. The integers ctj = 1,2, . . . , N r denote the 
replica indices. The integers m, n correspond to the dis- 
crete set of the Matsubara frequencies u> n — nT(2n + 1). 
The integers <jj = ±1 and tj = ±1 are spin and valley 
indices, respectively. The effective action is 



is invariant under the global rotations Qnm^rl ( r ) 

t^QferlWI" -1 ]^ in the combined spin-valley 
space for u S SU(4). The presence of the parallel mag- 
netic field and the valley splitting generates the symme- 
try breaking terms: 



S sb = iz s A s Tra z 



S vb = iz v A v Tr t z Q, (20) 



S = S a + S F + S s b + Svb + So, 



(16) 



where a z and r 2 are Pauli matrices in the spin and val- 
ley spaces, respectively. The Q-independent part of the 
action read o 5 ' 23 

So = -2nTz Trr ? A+^ J d 2 v ( X " S A 2 + X ° A 2 ) . (21) 

with x2 v — 2z SiV /tt being a bare value of the spin (valley) 
susceptibility. 



C. jT-algebra 

The action (Tl6|) involves the matrices which are for- 
mally defined in the infinite Matsubara frequency space. 
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In order to operate with them we have to introduce a 
cut-off for the Matsubara frequencies. Then, the set of 
rules which is called JF-algebra can be established.— At 
the end of all calculations one should tend the cut-off to 
infinity. 

The global rotations of Q with the matrix exp(i%) 
where x — J2a n Xn-^n pl & y the important role i 23 ' 24 For 
example, JF-algebra allows us to establish the following 
relations 

sp/;y*Qe-« = sp IZQ + 2in X °L n , 
tr V e^Qe-^ = tr V Q + £ in( X %Xt% sp I^Q^ 

cm 

- X> 2 (xX£(X-n)££S (22) 

an 

where sp stands for the trace over replica and the Mat- 
subara frequencies. 



III. SU(A) SYMMETRIC CASE 

A. jF-invariance 

At short length scales L <C L s , L v where L s<v — 
■\f Vxx/ (16z StV A StV ), the symmetry breaking terms S 3 b 
and S v b can be omitted and the effective theory becomes 
SU(4) invariant in the combined spin- valley space. Then, 
Eqs. (jTTJ) and (fT5)) should be supplemented by the impor- 
tant constraint that the combination z + T2 — 4F remains 
constant in the course of the RG flow. Physically, it cor- 
responds to the conservation of the particle number in 
the system^ In the special case of the Coulomb or other 
long-ranged interactions which are of the main interest 
for us in the paper the relation 

z + r 2 - 4r = (26) 



D. Physical observables 

The most significant physical quantities in the theory 
containing information on its low-energy dynamics are 
physical observables a' xx , z', and z' s v associated with the 
mean-field parameters a xx , z, and z SlV of the action (|16p . 
The observable a' xx is the DC conductivity as one can ob- 
tain from the linear response to an electromagnetic field. 
The observable z' is related with the specific heati^I The 
observables z' s and z v determine the static spin (x' s ) and 
valley (x' v ) susceptibilities of the 2D electron system£>2£ 
as x's v = 2z' SA] /ir. Extremely important to remind that 
the observable parameters a' xx , z' s v and z' are precisely 
the same as those determined by the background field 
procedure^ 

The conductivity a' xx is obtained from 



holds. With the help of Eqs. (|22[) . one can check that 
Eq. ([26]) guarantees the so-called ^"-invariance^, of the 
action S a + Sf under the global rotation of the matrix 



Q(r)e" 



(27) 



Here, Xn * s the unit matrix in the spin-valley space. In 
virtue of Eq. (|2"6"|) , it is convenient to introduce the triplet 
interaction parameter 7 = V2/ z such that T = (1+7)2/4. 
We notice that the triplet interaction parameter is related 
with Fq as 7 = -F$/(l + Fq ) . 



<4,(*«>») = -^( tr K= Q\[I-nM) 



(23) 



64Dn 



dr'((tr/£Q(r)VQ(r) tr J« n <3(r')VQ(r'))) 



after the analytic continuation to the real frequencies, 
ioj n — ► u> + i0 + at u> — > 0. Here, D = 2 stands for the 
space dimension, and the expectations are defined with 
respect to the theory lfT6|) . 

A natural definition of z' is obtained 2 ^ through the 
derivative of the thermodynamic potential £1 per the unit 
volume with respect to T, 



1 



d n 



2-7T tr r] A dT T ' 
The observables z' s v are given as 

7r d 2 Q 



B. Perturbative expansions 



To define the theory for the perturbative expansions 
we use the "square-root" parameterization 



(28) 



The action (11611 can be written as the infinite series in 



(24) the independent fields w%%:%% and w^gggg*. We 
use the convention that the Matsubara frequency indices 
with odd subscripts ni,ri3, . . . run over non- negative in- 
tegers whereas those with even subscripts n-2, rii, . . . run 
over negative integers. The propagators can be written 

(25) in the following form 
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a 1 a 2 -,<7 1 ,a 2 ( \ fa 4 a 3 -,a 4 a 3 ( \\ _ m ia3 J«2a4J ) X<?l<?3 X<72<?4 X X X D(u]^) 



32nTzj 



x£> p (o;i2)-D*(a;i2) 



(29) 



where W12 = w ni — uj n2 and 

^2; a; 



(30) 



[DIM]- 1 = p 2 



ig(z + r 2 )w„ 



C. Relation of z SjV with z and 7 

The dynamical spin susceptibility Xsi^, p) can be ob- 
tained froroi 

Xs(iv n , P) = X° s - Tz 2 s {tr iy z Q(p) tr /f „^Q(-p)> 

(31) 

by the analytic continuation to the real frequencies, 
ia; n — > lo + i0 + . Similar expression is valid for the val- 
ley susceptibility. Evaluating Eq. (f3Tj) in the tree level 
approximation with the help of Eqs. ()29|1 . we obtain 



,. , 2z s / 16z s w„ t 
X s («w„,p) = 1 D Juj n ) 

7T \ OV™ 



(32) 



In the case A s — A v — the total spin conserves, i.e., 
x(^,p = 0) = 0. In order to be consistent with this 
physical requirement, the relation 



z s = z + r 2 = z{\ + 7) 



(33) 



should hold. Similarly, the total valley isospin conserva- 
tion guarantees that 



z + T 2 = 2(1 + 7). 



(34) 



Being related with the conservation laws, Eqs. (|33|) and 
(|3"4")) are valid also for the observables: 



z' s = z' v = z'{l + i). 



(35) 



Therefore, three physical observables <j' xx , z' and 7' com- 
pletely determines the renormalization of the theory (TT 
at short length scales L <c L s , L v . 



D. One loop renormalization group equations 

As is shown in Ref . Q , the standard one-loop analysis 
for the action S a + Sf yields the following renormaliza- 
tion group functions that determine the zero-temperature 



behavior of the observable parameters with changing of 
the length scale L 



(36) 
(37) 
(38) 



d<Txx 


= /3a = 


- 2 -H + 

TT 


d;~) 


= /? 7 = 


(1+7) 2 
KOxx 


d\n z 


= lz = 


157-1 




KOxx 



Here, 7(7) = 1 - (1 + 7~ 1 )ln(l + 7), f = lni/Z and 
we omit prime signs for a brevity. Physically, the mi- 
croscopic length £ is the mean-free path length. It is the 
length at which the bare parameters of the action (Til)]) 
are defined. Renormalization group Eqs. (|36]) - ([38)) are 
valid at short length scales L <C L s , L v . 

As is well-known^ solution of the RG Eqs. |[3"6 > ([3T] ) 
yields the dependence of the resistivity p = 1/ ir<j X x on £ 
which has the maximum point and 7(£) dependence that 
monotonically increases with £. 



IV. St/(2) x 5(7(2) SYMMETRY CASE 

A. Effective action 

In this and next sections we assume that the spin split- 
ting is much larger than the valley splitting, A s ^> A„. 
Then, at intermediate length scales L s <c L -C L v the 
symmetry breaking term S s b becomes important. In the 
quadratic approximation it reads 



dr 



n 3 ,Tj 



(C2 



1/ nin2;riT2 H2ni;T2Ti 



(39) 

Hence, the modes in Q"^^ 1 ^ with o\ ^ <y% acquire a 
finite mass of the order of z s A s and, therefore, are neg- 
ligible at length scales L 3> L s . As a result, Q becomes 
a diagonal matrix in the spin space. Then, the spin sus- 
ceptibility has no renormalization on these length scales, 
i.e, 



dz s 
dt 



0, L S <^L<_L V . 
Let us denote Q***** 1 



(40) 



[QSL nn ]±- Then, the ac- 
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tion (fTBl) becomes S = S a + Sf + S v b where 



(41) 



and 



S F = AnTzY^ ( d 2 vtvrj{Q a — A) (42) 
+ nTfd 2 rJ2 ]T r aia2 trI%Q ai txr n Q a2 

an fTi,(72=i 

- nTT 2 J d 2 rJ2 ^ (tr J«Q CT ) ® (tr /^Q,). 

an cr— ± 

Now, the symbol tr stands for the trace over replica, the 
Matsubara frequencies, and the valley indices whereas 
Tr = J ePrtr. The action ((42|) corresponds to the fol- 
lowing low energy part of the Hamiltonian describing 
electron-electron interactions: 



Hint = r / (,V 



[J2 p ai T^p^ 



m a r t m a 



(43) 



It is worthwhile mentioning that Eq. (|4"3")l is in agreement 
with the ideas of Ref. 001 . 

The symmetry breaking part reads 



S sb = iz v A v ^2 / d 2 vtvT z Q 



(44) 



At length scales L ^ L s , the couplings T (T1<72 are all equal 
to each other, T crirT2 (L ~ L s ) = T = (z+T 2 )/4. However, 
the symmetry allows the following matrix structure of T: 



r | r++ r + _ 
1 1 r +H 



(45) 



As we shall see below, this matrix structure is consis- 
tent with the renormalization group. Physically, T++ 

and describe interactions between electrons with the 

same and opposite spins, respectively. 

The action (|4"Tj) and ([42]) is invariant under the global 



rotations [Q^ ;nT J CT (r) -> ^^[Q^rJ^r)^" 1 ]^ 
in the valley space for u a G SU(2). In order to preserve 
the invariance under the global rotations 

Q±(r) ^ e l *g±(r) e - 1 *, * = ( 46 ) 



where is the unit matrix in the valley space, the fol- 
lowing relation has to be fulfilled 



z + T 2 - 2r , 



2I\ 



(47) 



Physically, this equation corresponds to the particle num- 
ber conservation and is completely analogous to Eq. (|2"6")l . 



B. Perturbative expansions 

In order to resolve the constraint Q\ — 1 we use the 
"square-root" parameterization: 



Q± = W ± + KJl-Wi 



(48) 



Then, the action (|4T|) and (f4"2"|) determines the propaga- 
tors as follows 



where 



32 * 

HrtM = — 2><r«r>, (49) 



^ 1 ,n 3 ^ T1T3 <5 T2T4 £> g (a;i2,T 1 ,r 2 ) 



32ttT 



32ttT 



f5 Q1Q2 5 T1T2 ,5 T3T4 ^(wi 2 )D^i 2 ) 



(50) 



with 
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K(w n )]- 1 = g 2 + (z + r 2 -2I>„ (51) 



D-\u n , 71,7a) = D" 1 ^) + *^^(r 1 - r 2 ), (52) 

[Dj(w n ,Ti,7a)] -1 = [D*K)]- X + i— (n - r 2 ). 

(53) 

In the same way as in Sec. IIII CI the conservation of the 
total valley isospin guarantees the relation z v — z + T 2 . 
The conservation of the z-component of the total spin, 
p + — p~ , implies that z s = 4T_| (see Eq. (|3~TjO . There- 
fore, 



dhr 4 







d£ 

for the length scales L s <C L <C £ t 



(54) 



C. One-loop approximation 

Evaluation of the conductivity according to Eq. ([23 
in the one-loop approximation yields 



v'xxiiun) = cr xx + 



2 8 7T 



P 

P uj m >0 



1 DjK) 



',D p (LO m + LU n ) 

[J2 (TD s p (u; m ))^-4r 2 D p (u 



Hence, we find 



<T=± 



<T'xx( iu >n) = Oxx + [ P 2 T V min{ — , l}D p (u m ) 



(55) 



P Wm>0 



xD p (u m + uo n ) zD s p {uo m ) - QT 2 D t p {u) m ) 

-(z + 2r 2 -4r ++ )^*K l )l (56) 
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where 



G4 



-uj„T+-. 



(57) 



Performing the analytic continuation to the real frequen- 
cies, iu> n — ► u> + i0 + in Eq. (|5rj|) , one obtains the DC 
conductivity in the one-loop approximation: 

2 8 7r /" f 00 

P 2 dwDl{w) zD s p (u) 



-6r 2 D^( w ) - (z + 2T 2 - 4r ++ )/3» 



(58) 



In order to compute z' and z(, we have to evaluate the 
thermodynamic potential Q in the presence of the finite 
valley splitting A„. In the one- loop approximation we 
find 

,dtt/T 



T 2 ^1LL = 8N r T Y 
dT ^ 



UJ n >0 



+ — f f2r+_^*(o; n ) 

Vxx Jp L 

{z + r 2 )Dl (w n ) + (z + r 2 ) fpK,n,T 2 ) 

£>p(a; n ,Ti,r 2 )ll. (59) 



Following definitions (|24p and (|25p of the physical ob- 
servables, we obtain from Eq. (fS"9"|) 



— (2r 2 -r ++ ) / d p (o) 

®XX J r> 



(60) 



and 



1 + 47T 
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(z + T 2 )T / f^>„) 



O)„>0 



(z + Ts)^ 3 ^) 



(61) 



We mention that the results (|58|) , (|60| , (fBTj) can be ob- 
tained with the help of the background field procedure^ 
applied to the action (|4"Tj) - (|42|) . 



D. One loop RG equations 

Using the standard method,— we derive from 
Eqs. (|58|) . (|60| and foTj) one-loop results for the RG equa- 
tions which determine the T — behavior of the physical 
observables with changing the length scale L. It is conve- 
nient to define j v — T 2 /z and j s = —1 + 4r_| j z. Then, 

for D = 2 we obtain 



d<Jxx 


-i [1 + 6/(7*) 4 

7T 


■/(7.)] 


(62) 




(1 + 27„ - 

ncr xx 




(63) 


da ~ 


(1 - 67„ - 

TCTxx 


-7.) 


(64) 


dlnz 
d£ 


— [1 - 6 7 „ 

K(Jxx 


-Is]- 


(65) 




FIG. 1: The projection of the RG flow in the three dimen- 
sional parameter space (a xx , j v , j 3 ) onto (7*,, 7s) plane for the 
5(7(2) x 517(2) symmetry case (Eqs. ([62 [l -[[6l } '). Dots denote 
the line at which 1 + 6/(71,) + /(7s) = 0. The dashed line 
indicates the line 7„ = 7 S (see text). 




FIG. 2: Schematic dependence of the resistance p — l/(na xx ) 
on £. Curves a, b, and c corresponds to the flow lines a, b, 
and c in Fig. [T](see text). 



Eqs. ((62 |) - (f65|) constitute one of the main results of the 
present paper and describe the system at the intermedi- 
ate length scales L s C L « i„. We mention that the 
length scale I involved in £ = In L/l is now of the order 
of L s . 

In Figure[T]we present the projection of the RG flow in 
the three dimensional parameter space (o'xxtlfvi'Ys) onto 
ilvi Is) plane. There is the unstable fixed point at 7^ = 
and 7s = 1. However, for the physical system considered 
the fixed point is inaccessible since an initial point of the 
RG flow is always situated near the line j v = j s . As 
shown in Fig. [21 there are possible three distinct types of 
the p(£) behavior for such initial points. Along the RG 
flow line a (Fig. [T]) that crosses the curve d described by 
the equation l + 6/(7„) + /(7 s ) = the resistance demon- 
strates the metallic behavior: p decreases as £ grows. If 
we move along the RG flow line b which intersects the 
curve d twice, then the resistance develops the minimum 
and the maximum. At last, the resistance on the RG flow 
line c which has single crossing with the curve d has the 
maximum. Remarkably, in all three cases, the behavior of 
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the resistance is of the metallic type for relatively large L. 
The reason of this metallic behavior can be understood 
from the following arguments. At large £, the coupling 
7„ flows to large positive values whereas "f s — > — 1. Then, 
r+_/r ++ - l/7„ < 1 and the RG Eqs. (|62] ) -(|65 |) trans- 
forms into equations for the single valley system with the 
conductance equal a xx /2. The metallic behavior of this 
system is well-known^ 



V. COMPLETELY SYMMETRY BROKEN CASE 

A. Effective action 

At the long length scales L^> L v the symmetry break- 
ing term becomes important. In the quadratic ap- 
proximation it reads 



S v b = 



tz v A v 



Oij , (7 

rij ,Tj 



( T 2 T l) [ w nln2;T 1 T 2 \v[ W n2nl-,T 2 T 1 ]'7 ■ 



(66) 

Hence, the modes in [Q"m TlT2 ]cr with t\ ^ r 2 acquire 
a finite mass of the order of z v A v . Therefore, they are 
negligible at long length scales L 3> L v . As the result, the 
matrix Q becomes diagonal matrix in the valley isospin 
space. The valley susceptibility remains constant under 
the action of the renormalization group on these length 
scales: 



dz v 

— - = 



L > Ly. 



Let us define 



Qf = {[Q a 11 p } + AQ-1 

Then the action S = S a + Sf reads 



~32~^ 



reads 
J d 2 rtr(VQj) 2 



and 



S F = nT I d'rJ2J2 tTl nQ^trIl n Q k 

j.k an 



d rtrr/Qj, 



(67) 
(68) 
(69) 

(70) 
(71) 



where 



r = 



/r++ - r 2 


r++ 


r+- 


r+- 


r++ 


r++ - r 2 


r+- 


r+- 


r+- 


r+- 


r++ - r 2 


r++ 


V r + _ 


r+- 


r++ 


r++- 



r +- \ 



r 2 / 



(72) 

Initially, at the length scale of the order of L v , the cou- 
pling r ++ = r ++ and r_| = r_| However, more gen- 
eral structure (|72p is consistent with the renormalization 



group. It is worthwhile to mention that if the matrix 
r is diagonal then the theory (f6"9"|) and (|71[) would in- 
clude four copies of the singlet U(l) theory studied in 
Refs. [7.8]. The action (fTTj) corresponds to the follow- 
ing low energy part of the electron-electron interaction 
Hamiltonian: 



gg'\tt ! 

P% = (73) 



In order to have the invariance under the global rota- 
tions 



Qj 



x 



E aja 



(74) 



(75) 



the following relation has to be fulfilled 

z + r 2 -r ++ -f ++ = r+_ + f + _. 

B. Perturbative expansions 



As above, in order to resolve the constraints Q? = 1, 
we shall use the "square-root" parameterization for each 



Qj'- Qj = w j + A y 1 - Wf. Then, the propagators are 
defined by the theory ([69|) and (fTTj) as 

([■"Cnl il)]j \ w ZT/' TiTS (-Q)]k) = —V jk , (76) 



V = 6 aia3 5 a2a4 5„ 



6 nx , n3 D q (u 12 ) + ^lt6 a ^ 



y.D q {Lo 12 )b c q {Lo 12 ) 



where 
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KK)]" 1 =q 2 + {z-T)u n 



(77) 



(78) 



The conservation of the ^-components of the total 
spin, J2ctt' j P'tj an d the total valley isospin, ^2 aT Tp%, 
implies (see Eq. ([31])) that z s = 2r + _ + 2f + _ and 

z v = 2T ++ + 21^ Since, for L ^> L v both z s and 

z v are not renormalized, we obtain 



dT, 



di 



= 0. 



(79) 



Since, both 1% and r_| coincides at the length scales 

L ~ L v and they are not renormalized we shall not dis- 
tinguish and r_| from here onwards. If we intro- 
duce "f s and 7n such that r_| = r_| = z(l + 7 s )/4 and 

r ++ = z(l + 27„ — 7s)/4 then both 7 S and j v coincide 
with the corresponding couplings of the previous sections 
at the length scales L ~ L v . 
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C. One-loop approximation 

Evaluating the conductivity with the help of Eq. ([231 
in the one-loop approximation, we find 



a xx( iuJ n) = Oxx + 



— / p 1 > mm < , 1 > 

DaxxJ/ J-^ J 



D p (uj m + uo n )D p (uj m ) ^ (r-Dp(^m) 



(80) 



Hence, 



2 8 7 



mm <, , 1 



'P w m >0 

D p {u) m + Lu n )D p (uj m ) Dp(uj m ) - 2j v D p (uj m ) - j s D p (u m ) 

(81) 



where 



32 



[5*K)]- 1 = g 2 + — (r ++ + f 



(82) 




FIG. 3: The projection of the RG flow in the three dimen- 
sional parameter space (a xx , 7„, 7 S ) onto (7^, 7s) plane for the 
completely symmetry broken case (Eqs. (|86|) - (|88| )'l. Dots de- 
note the line at which 1 + 2f(-y v ) + /(7s) = (see text). 



Performing the analytic continuation to the real frequen- 
cies in Eq. (j8Tj) . we find 



2 8 7TZ 



a xx ~ a xx 



p 2 / dLoDliu) D s (lo) 



'xx J p JO 
V 



(83) 



As in the previous Section, in order to compute z' we 
evaluate the thermodynamic potential in the one-loop 
approximation. The result is 



dT 



W„>0 



1 



Hence, we obtain 

16 



x =z + — (r 2 

®XX 



D p (0). 



(84) 



(85) 



We mention^ that the results (JTHD, (JHSJ), and 
can be obtained with the help of the background field 
procedure applied to the action ([6"9"]) - ([7T]) . 



D. One loop RG equations 



Equations ([79]), ([83]) and ([85]) allow us to derive the 
following one-loop results for the renormalization group 
functions which determine the T = behavior of the 
physical observables with changing the length scale L 



FIG. 4: Schematic dependence of the resistance p — l/(na x 
on £ along the flow line a in Fig. [3] (see text). 



(D = 2): 



= -![! + 2/(7,) +/( 7 .)] 

-7T = (1 - 27^ ~ 7s) 

at, ttct xx 

^ = ^(1-^-7.) 

CI? IT (Jr. 



din i 



'XX 

1 

ira x 



[1 - 27„ - 7 S 



(86) 
(87) 
(88) 
(89) 



The renormalization group equations ([86 ]) - ([89] ) constitute 
one of the main results of the present paper. We mention 
that the length scale I involved in £ = InL/l is now of 
the order of L v and Eqs. ([8rj]) - ([89"]) describe the system 
at the long length scales L 3> L v . 

The projection of the RG flow for Eqs. po]) -(|88 ]l on 
the j v - 7 S plane is shown in Fig. [3] There exits the 
line of the fixed points that is described by the equation 
27u + 7s = 1 • If the initial point has large j v or 7 S then 
the RG flow line crosses the curve that is determined by 
the condition 1 + 2/(7„) + /(7 S ) = 0. Therefore, the p(£) 
dependence along the RG flow line develops the minimum 
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A T<« 




FIG. 5: The schematic p(T) dependence in the case of zero 
parallel magnetic field. See text 



FIG. 6: The schematic p(T) dependence in the presence of 
both spin and valley splitting in the case A s < A v . For the 
opposite case, the behavior will be similar. See text 



and will be of the insulating type as is shown in Fig. 3J 



VI. DISCUSSIONS AND CONCLUSIONS 

The renormalization group equations discussed above 
describe the T — behavior of the observable parameters 
with changing of the length scale L. At finite tempera- 
tures T > <r xx /(zL 

tampie) wu ^re L S ampie is the sam- 
ple size, the temperature behavior of the physical ob- 
servables can be found from the RG equations stopped 
at the inelastic length L; n rather than at the sample 
size. Formally, it means that one should substitute 
£t = ^ In a xx /(zTl 2 ) for £ in the RG equations with £y 
obeying the following equation 31 



1 



1 din z 

2 d£ ' 



(90) 



Having in mind Eq. (|90| . we find that the T-behavior 
of the resistivity at B = is described by Eqs. ([55)) and 
([57]) for T > A„ and Eqs. ([52 ]) -([Si ]) with interchanged 
7„ and 7 S for T <C A„. In what follows, we assume 
that A v < Tm2x where T^/ax denotes the temperature of 
the maximum point that appears in p{T) according to 
the RG Eqs. ([55)1 and ([57)1 . Our assumption is consistent 
with the experimental data in Si-MOSFET where, for ex- 
ample^ the valley splitting is of the order of hundreds 
of mK and T^lx is about several Kelvins 3 - 3 -. Then, de- 
pending on the initial conditions at T ~ 1/t two types of 
the p{T) behavior are possible as is shown in Fig.0 The 
curve a represents the typical p(T) dependence that was 
observed in transport experiments on two- valley 2D elec- 
tron systems in Si-MOS samples 2 and n-AlAs quantum 
welli^ Surprisingly, the other behavior with the two max- 
imum points is possible, as illustrated by curve b in Fig. [5] 
So far, this interesting non-monotonic p(T) dependence 
has been neither observed experimentally nor predicted 
theoretically. At very low temperatures T <C A„, the 
metallic behavior of p(T) wins even in the presence of 
the valley splitting. 



In the presence of the sufficiently low parallel magnetic 
field A s < Tmax, the p(T) behavior of three distinct types 
is possible as plotted in Fig. [(51 In all three cases, the 
p{T) dependence has the maximum point at tempera- 
ture T = Tn2x and is of the insulating type as T — > 0. 
As follows from Fig. O in the intermediate temperature 
range, when T is between A s and A„, the metallic (curve 
a) , insulating (curve b) and nonmonotonic (curve c) types 
of the p{T) behavior emerge. As a result, there has to 
exist the p(T) dependence with two maximum points in 
the presence of B. 

For high magnetic fields such that A 



> TgL the 



maximum point at T = Xll/ax is absent, and two types 
of the p(T) behavior are possible as is shown in Fig. [7J 

If TmaL < A„, then the dependence of the resistivity 
is monotonic and insulating, see the curve a in Fig. [7J 
Here, T^Ix denotes the temperature of the maximum 
point that appears in the resistivity in accord with the 
RG Eqs. ([12]) and ([55)1. In the opposite case t£2 > A„, 
a typical p(T) dependence is illustrated by the curve b 
in Fig. [7J Therefore, if the valley splitting is sufficiently 
large, i.e., A^ > T^Il, then the monotonic insulating 
behavior of the resistivity appears in the parallel mag- 
netic field which corresponds to A s ~ Tm2x. This is the 
case for the experiments on the magnetotransport in Si- 
MOSFET J^ 3 . However, if the valley splitting is small, 
An < Tm&x, then the maximum point of the p(T) de- 
pendence survives even in high magnetic fields but shifts 
down to lower temperatures. 

In addition, to interesting T-dependences of the resis- 
tivity, the theory predicts strong renormalization of the 
electron-electron interaction with temperature. In order 
to characterize this renormalization, we consider the ra- 
tio Xv/Xs of valley and spin susceptibilities. In Figure [51 
we present the schematic dependence of Xv/Xs on T for a 
fixed valley splitting but with varying spin splitting. At 
high temperatures, T > A„, A s the ratio of the suscep- 
tibilities equals unity, Xv/Xs = 1. At low temperatures, 
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p 




a 



the temperature behavior of such physical observables as 
the resistivity, spin and valley susceptibilities in 2D elec- 
tron liquid with two valleys in the MIT vicinity and in 
the presence of both the parallel magnetic field and the 
valley splitting. First, we found that the metallic be- 
havior of the resistivity at low temperatures survives in 
the presence of only the parallel magnetic field or the 
valley splitting. If both the spin splitting and the val- 
ley splitting exist then the metallic p-dependence crosses 
over to insulating one at low temperatures. Second, we 
have predicted the existence of the novel, nonmonotonic 
dependence of resistivity at zero and finite magnetic field 
in which the p(T) has two maximum points. It would be 
an experimental challenge to identify this novel regime. 
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FIG. 8: The schematic dependence of the ration Xv/Xs on 
temperature: a) for A s < A„; b) for A s = A„; c) for A s > 
A v . The temperature scales T a ,c = A a . 



T <C A. u , A s , we find 



Xs 



< 1 , A s < A v 
= < 1 , A s = A„ 
> 1 , A s > A^. 



(91) 



Therefore, the ratio Xv/Xs at T — * is sensitive to the 
ratio A„/A s . This can be used for the experimental de- 
termination of the valley splitting in the 2D electron sys- 
tem. 

Finally, we remind that we do not consider above 
the contribution to the one-loop RG equations from the 
particle-particle (Cooper) channel. It can be shown (see 
Appendix) that neither the spin splitting nor the val- 
ley splitting does not change the "cooperon" contribu- 
tion to the RG equations in the one-loop approximation. 
Therefore, the Cooper-channel contribution to the RG 
equations discussed above can be taken into account by 
the substitution of 1 + 2 for 1 in the square brackets of 
Eqs. {35]), ||5JJ) and §^ . The Cooper-channel contribu- 
tion does not change qualitative behavior of the resis- 
tivity, and the valley and spin susceptibilities discussed 
above. 

To summarize, we have obtained the novel results on 



APPENDIX A: "COOPERON" CONTRIBUTION 
TO THE CONDUCTANCE 

We start from the standard equation for the 
"cooperon" 



1 



where the "impurity ladder" is given as 



(Al) 



reri ,(T2 ;tiT2 ( „\ _ 

CT30"4;T3T4 Vi) 



(~iR 

2-KVT, I ^W>™ 



(P+)G* 3 a 4 ;T 3 T 4 (P-) 
(A2) 



Here, p± = p ± q/2 and the impurity averaged Green 
functions are given as 



p- 



2m 

+A^ sgriTi ±i/2r, 



fi + A s sgncri 
(A3) 



Performing integration, we find for q — ► 



7-CTl ,o- 2 ;tit 2 

J (T3 (T4;T3T4 



l - Dq 2 T t 



-iA s Ti(sgna 3 - sgncri) - iA v Ti(sgnT 3 - sgnn) 



(A4) 
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Next, solving Eq. (|Aip . we obtain 



we find 



Gil 



1,<T2;T1T2 ( \ JTCTlfJ2 rCT3(T4 TT1T2 rT3T4 



2-KVTf 



Dq 2 (A5) 



-iA s (sgn cr 3 - sgncri) + iA„(sgnr 3 - sgnri) 



The interference correction to the conductance is given 
as 



Sa xx = -- f G R (-p) 

n J P 

xC R (ri\C A 

Using the result: 



f 3 CT 4 ;T3 T4 G(Ji 0-5 JT4T5 

(-p) (A6) 



(P) 



d 2 g 



/-^CT2 CT 3 ;T2T3 

2 cr 5 o-6 ;t 5 t 6 



(?)• 



(p)G (p) 



Sa xx = -J2 6*^5™ [l + A 2 T 2 (sgn<7i _ sgn<72 )2 

d 2 q 



+A t ,r l (sgnn - sgnr 2 ) 



D 



(27T) 



TV 



-zA s (sgncri - sgn cr 2 ) + iA„(sgnri - sgnr 2 ) 



9 ' 



(A8) 



47rz/T? J* 73 CT5 <5 g2 J6 <5 T3 Tj (? T2 Tb Therefore, neither the spin splitting nor the valley split- 

1 + A2r 4 2 (sgn<T2 — sgn<73) 2 + A 2 r 2 (sgnr2 — sgnrs) 2 ' ting affect the Cooper channel (interference) contribution 

(A7) to the conductance. 
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